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Abstract 

We study the stability of the Einstein static universe in the Hof ava-Lifshitz (HL) gravity and 
a generahzed version of it formulated by Sotiriou, Visser and Weifurtner. We find that, for 
the HL cosmology, there exists a stable Einstein static state if the cosmological constant A is 
negative. The universe can stay at this stable state eternally and thus the big bang singularity 
can be avoided. However, in this case, the Universe can not exit to an inflationary era. For 
the Sotiriou, Visser and Weifurtner HL cosmology, if the cosmic scale factor satisfies certain 
conditions initially, the Universe can stay at the stable state past eternally and may undergo 
a series of infinite, nonsingular oscillations. Once the parameter of the equation of state w 
approaches a critical value, the stable critical point coincides with the unstable one, and the 
Universe enters an inflationary era. Therefore, the big bang singularity can be avoided and a 
subsequent inflation can occur naturally. 
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I. INTRODUCTION 



In the standard cosmological model, the existence of a big bang singularity in the early 
universe is still an open problem. In order to resolve this problem, Ellis et al. proposed, 
in the context of general relativity, a scenario, called an emergent universe l|, l2|. In 
this scenario, the space curvature is positive and the Universe stays, past eternally, in 
an Einstein static state and then evolves to a subsequent inflationary phase. So, the 
Universe originates from an Einstein static state, rather than from a big bang singularity. 
It is also worth noting that an Einstein static state as the initial state of the Universe is 
also favored b. the eut.opy consideration, g However, the Einstein statie universe ,n 
the classical general relativity is unstable, which means that it is extremely difficult for 
the Universe to remain in such an initial static state in a long time due to the existence of 
perturbations, such as the quantum fluctuations. Therefore, the original emergent model 
does not seem to resolve the big bang singularity problem successfully as expected. 

Since in the early epoch, the Universe is presumably under extreme physical condi- 
tions, new effects, such as those stemming from quantization of gravity, or a modification 
of general relativity or even other new physics, may become important. As a result, the 
stability of the Einstein static state has been examined in various cases 4l-ll6l]. For exam- 
ple, the emergent scenario within the frameworks of loop quantum gravity and braneworld 



cosmology have been discussed in Refs. 



where it was found that a successful model 



can be obtained, while the stability in the presence of vacuum energy corresponding to 
conformally invariant fields has been studied and a nonsingular emergent cosmological 
model was reconstructed j4|. In f{R) gravity, it was found that the Einstein static state 
is stable under homogeneous perturbations [9|, but this stability is broken by adding in- 
homogeneous perturbations 10|. In f{G) gravity, the stability of the Einstein static state 



agan.t homogeneous perturbations has bjn analyzed in Ref. where G rs the Guass- 

Bonnet term. In addition. Barrow et al. [12] found, with the covariant techniques, that the 
Einstein static state is stable for small inhomogeneous vector and tensor perturbations, 
as well as for adiabatic scalar density inhomogeneities with > 0.2. 

Recently, motivated by the Lifshitz theory in solid state physics, Hofava proposed a 
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power-counting renormalizable theory of gravity, called Hofava-Lifshitz (HL) gravity 17 1. 
In the ultraviolet (UV) limit, HL has a Lifshitz-like anisotropic scaling between space and 
time characterized by the dynamical critical exponent z = 3 and thus breaks the Lorentz 
invariance, while in the infrared (IR), it flows to 2; = 1. So, it is expected to reduce to the 
classical general relativity gravity theory in the low energy limit. Applying the HL gravity 
to cosmology, it has been found, in a nonflat universe, that the Friedmann equation is 
modified by a ^ term 18|-|20[|. where a is the scale factor. The cosmological perturbations 
with the HL gravity were studied in Refs. 



21 



and the results showed that a scale 



invariant superhorizon curvature perturbation could be produced without inflation. In the 
original HL gravity, Hofava assumed two conditions: detailed balance and projectability. 
More recently, Sotiriou, Visser and Weifurtner (SVW) [31] suggested to build a general 
HL theory with projectability but without detailed-balance conditions. For a spatially 
curved Friedmann-Robertson- Walker universe, the SVW generalization gives an extra 
correction term and modifies the coefficient of the ^ term in the Friedmann equation as 
compared to the HL theory. Let us note here that some other issues in HL gravity have 



been dealt with in Refs. 



32 



35|. 



In this paper, we will discuss the stability of the Einstein static universe in the contexts 
of HL gravity and SVW HL theory, respectively. In the following section, we briefly review 
the HL and SVW HL cosmology. In section HI, we analyze the Einstein static solutions 
and discuss the stability of these solutions. Finally, in section IV, we present our main 
conclusions. 



II. THE HORAVA-LIFSHITZ COSMOLOGY 



In HL gravity, it is convenient to use the Arnowitt-Deser-Misner decomposition of the 
metric 

ds^ = -N^dt^ + gij{dx' + N'dt){dx^ + N^dt) , (1) 

where Qij is the 3-dimensional spatial metric, is the lapse function, A^* is the shift 
vector, and the coordinates scale as t — )■ i^t, — )■ ix\ In this paper, we only consider 
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the z = 3 case. Next, we first turn our attention to the implications of HL gravity in 
cosmology. 



A. The HL cosmology 

The action of HL gravity consists of kinetic and potential terms. The former is given 

by 

^fe = ^ y dtd^x^N{K,,K'^ - \K^) , (2) 

where Kij = ^{(jij—ViNj—VjNi) is the extrinsic curvature, K = g^^Kij, K^^ = g^'^g^'-Kki, 
and A is a dimensionless parameter. When A = 1, one recovers the kinetic part of the 
4-dimensional Einstein-Hilbert action. With the detailed-balance condition, the potential 
term has the form 

Sv= f dtd^x^N{(3CijC'^ + -f^RaVjRi + (RijR'^ + r]R^ + ^R + a) . (3) 

Here Rij is the three-dimensional spatial curvature tensor, R = g^^Rij, e^^^ is the anti- 
symmetric tensor with e^^^ = 1 and C*-' = ^Vk{Rj — ^^iR) is the Cotton tensor. The 
constants (3, 7, (, r], ^ and a are defined, respectively, as 

_ _kV^(1-4A) 
^~~2^' ^~2^' ^~ 8~' 32(1 -3A) ' 

/^V ^ 3kV ^2 

^ 8(1 -3A) ' 8(1 -3A) ' ^ ' 

where A is the cosmological constant, and fi and u are two coupling constants. In this 
case, the emergent speed of light becomes 

So, for the case where 3A — 1 > 0, a negative A is required in order to guarantee that 
the speed of light is real. Let us note that a positive A can be obtained by making an 



analytical continuation for parameters fi and u"^ hj ^ ifi and — i- — iw^ 18|. In 



addition, it was found in Ref. 



36| that a negative cosmological constant may disappear 



in the different geometries, plane symmetric spacetimes, for example. 



For a homogeneous and isotropic universe described by the metric 



ds^ = -dr + a\t) 



dr^ 



+ r^d^n 



1 — kr^ 

the Friedmann equation in HL gravity can be expressed as 

— [3X-1)H =p-a 



(6) 



K 



a" 



(7) 



Here k = 0, ±1, and p is the energy density of a perfect fluid in the universe, which 
satisfies the conservation equation 



p + 3Hp{l + w) = 



(8) 



where w = p/p is the equation of state. In the present papei\ a constant w is considered, 

fl Q 

which is a good approximation, since, as shown in Refs. a plateau potential is 

required in the past-asymptotic limit in the emergent scenario. It is easy to see that, for 
a spatially flat universe, this Friedmann equation is the same as that in general relativity. 
Now, we define two new constants 



a 



2 4 



(9) 



6(3A-1)' 48(3A- 1)(1 -3A) 

Apparently, (3 is positive if 3A — 1 > and negative if 3A — 1 < 0, whereas a is always 
negative if yU^ > 0. Using these newly defined constants, the above Friedmann equation 
for a closed universe can be reexpressed as 



H^ = /3p + 3aA^ ^ + ^ 



(10) 



which be further written as 

= /3p + 3aK^ 



2 1 

1- ^ + 



Aa2 A2a4 



Differentiating this equation with time and using the energy conservation equation, we 
have 



a 

a 



3w 



3(1 + 3w) 



Aa2 



3(1 -3w) 1 
2 



(12) 
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B. The SVW HL cosmology 



Sotiriou, Visser and Weifurtner generalized the original 
jectability but abandoning the detailed-balance conditions 
Friedmann equation has the following form 

3 



iL theory by keeping the pro- 



3l|. In this case, the modified 



^ 3 a2 a4 ^ a6 



(13) 



where A, Pi and /32 are coupling constants. Comparing the above equation with that in 
HL theory, one can see that the SVW generalization not only gives an extra correction 
term, but also modifies the coefficients of other terms. It is interesting to note that the 
co„ectio„ tenn ,„ay also .esu.t .o. a .adiat.o. fluid w,th,„ the UV .eg-^e 0. 
So, the influence of such a radiation fluid within the UV regime on the stability of an 
Einstein static state can be regarded as a specific case of the analysis to be carried out 
next. Defining two dimensionless constants, /3i = PiA and = ^2-^'^, we find that the 
modified Friedmann equation in SVW HL theory becomes 



1 , 2/3i ^ 4(32 



Aa2 A2a4 

Then, differentiating Eq. (IT^ with time, one has 



Ah 



(14) 



2(1--A 



3w)H' 



A 



, l + 3w 6w-2^ 12w 
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Aa2 



A2a4 



Ah 



B2 



(15) 



III. THE EINSTEIN STATIC SOLUTION 



The Einstein static solution is given by the conditions a = and a = 0, which imply 



a = aEs, H{aEs) = . 



(16) 
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A. The HL cosmology 



From Eq. (fT2|) . it is easy to that the Einstein static solution satisfies the following 
equation 

Solving this equation, one obtains two critical points 

Point A: i = A (18) 



and 



Po^nt B : ^ = -^ii±^A . (19) 
Substituting these critical points into Eq. (fTOl) reveals that Point A corresponds to 



and Point B to 



PA = (20) 



48a; .n 



If A is negative, Point A is physically meaningless since al;^ = x negative. The existence 



condition for Point B is —^73^ < 0, which leads to —1 < w < \. For a positive 



l-3w 



A obtained from an analytical continuation for parameters fi and u"^ hj fi ^ ifi and 
cj^ — )■ —ioj"^, which yields a positive a since /i^ becomes negative after the analytical 
continuation [refer to Eq. (Q], it seems that Point A exists and so does Point B if w 
satisfies the condition w < —1 or w > ^. However, in this case, the energy density 
corresponding to Point B, ps, becomes a negative since a is positive, which is meaningless. 
So, in the case of a positive A, only point A exists physically. 

In order to study the stability of these critical points, we introduce two variables 

Xi = a, X2 = a . (22) 

They obey the following equations 

xi = X2 , (23) 
7 



3w xl 9 .n,^ , o A /1 o \ 1 3a(l — 3w) 1 

^ + -aA^ 1 + w)xi - 3aA 1 + 3w) 

I Xl 2 Xl 2 x-{ 



Linearizing the system described by the above two equations near critical points, one can 
obtain the eigenvalues of the coefficient matrix, which determine the stability of these 
critical points. After some calculations, we get the eigenvalue 'i?^: 

Point A: = 12aA2 , (25) 

Point B : ^'^ = -12aA^ . (26) 

«L 

If "i?^ < 0, the corresponding equilibrium point is a center point, otherwise it is a saddle 
one. In order to analyze the stability of the critical points in detail, we now divide our 
discussions into two cases, i.e., A < and A > 0. 



1. A < 

In this case, a < and Point A is physically meaningless since a|;^ < 0. Therefore, we 
only discuss the stability of Point B, which exists under the condition — 1 < w < |. From 
Eq. (l26ll . it is easy to see that ■(9^ < since aA > 0. This means that point B is stable. 
In Fig. [H we plot the evolution of the scale factor with time and the phase diagram in 
space (a, a). This figure shows that the Universe can stay at the stable state eternally 
and may undergo a series of infinite, nonsingular oscillations about this point. Thus, the 
initial big bang singularity can be avoided. By numerical calculation, however, one finds 
that, when w is larger than | or less than —1, the Universe may undergo an accelerating 
expansion. It therefore appears that the universe may enter an infiationary phase from 
this stable point if the condition — 1 < w < | is violated. However, from Eq. ( fT9l) . we find 
that, once w evolves through —1 or |, the scale factor a becomes oo or 0. Therefore, in 
this case, the Universe is essentially stuck at the stable static state unless the scale factor 
becomes singular. 
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FIG. 1: The evolutionary curve of the scale factor with time (left) and the phase diagram in 
space (a, d) (right) for the case A < in Planck unit and with w = —0.90, A = —0.6, a = —1. 



2. A > 



A positive A can be obtained by making an analytical continuation of the parameters 
/i and by ^ ifi and u"^ — )■ iu'^ [l8]. The analytical continuation changes the sign of 
a and makes it a positive constant {a > 0), since it contains a /i^ factor. Now, critical 
point A can exist, but it is a saddle point since > 0. The critical point B is physically 
meaningless due to < as we have pointed out. Therefore, in this case, there is no 
stable Einstein static universe. 

A summary of the existence and stability of Points A and B is given in Table [H 



TABLE I: Summary of the critical points and their stability in the case of the HL cosmology 





A < 


A > 


Existence 


Stability 


Existence 


Stability 


Point A 


meaningless 






unstable 


Point B 


-1 <w < I 


stable 


meaningless 





B. The SVW HL cosmology 

For the SVW HL cosmology, we only consider the case of a positive cosmological 
constant (A > 0)^. From Eq. (fTSl) . one can see that the critical points are determined by 

^ The negative cosmological constant case can be treated similarly 
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the following cubic equation: 

, 1 + 3w 6w - 2 - 12w - 12 - 

When P2 = 0, the above equation simplifies to a quadratic one, which is similar with 
that in the HL theory, but not identical, since coefficients are different. Thus, now we 
separately discuss two cases, (32 = and /32 7^ 0. 

1. (32 = 

In this case, Eq. (|27|) reduces to 

N 1 + 3w Qw - 2 - ^ 

Solving this equation, one can obtain two critical points: 

1 _ 1 

~ 4(3w;- l)/3i 



Point C ^^ = — __[i + 3u;- -y/(l + 3w)2-8/3i(3w2 + 2w;-l)] , (29) 



Point D = — [1 + 3w + V (1 + 3w)2 - 8/3i (3^2 + 2w- 1)1 . (30) 

When /3i = |, critical points C and D reduce to = |A and ^ = ^^^(^A, respectively, 
which are the same as that in the HL cosmology [given in Eqs. (fT8|[T9l) ] after a redeffnition 
of the cosmological constant as |A. It follows, from Eqs. fl29l [30]) . that when 



^ (l+y)^ (31) 

8(3w-l)(w; + l)' ^ ' 

two critical points, C and D, coincide, and thus, in this case, there is only one critical 
point, 

i = 4(I;7TT)^ • 

The energy density p, at critical points C and D, is 
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where is given by (l29l) or (l30l) . In order to ensure these critical points to exist with 

"■Es 

a physical meaning, it is required that a^^ > and p{aEs) > 0. This yields a region of 
existence in the {w,/3i) parameter space, 
• For Point C: 



For Point D: 



w < < ^Si < ^ , 

1 ^ ^1 (l + 3w)2 

~3-^-3' 8(3^2 + 2w-l] 



1 1 {i + SwY 

— < w < -, , , 

3- 3 8(3ti;2 + 2ti; - 1) 



< /3i < . 



(34) 



(35) 



In Fig. [21 we show the regions of existence in the (w, /3i) parameter space for both critical 
points C and D. 



2 
1 
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-3 
-4 
-5 



1 


-1 
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-5 



FIG. 2: Regions of existence in the (t(;,/3i) parameter space. The left panel shows the existence 
region for Point D, while the right panel for Point C . 

With the same method as that used in the HL theory, we find the eigenvalue of critical 
points C and D, which can be expressed as 

1 + 3w 3(6ry - 2)/3i 



;i + ^) 



(36) 



There is no point in the existence region in the parameter (ti?, /3i) space for Point C which 
gives rise to a negative '{p' . Hence Point C is always unstable. For critical Point the 
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region of stability and existence is 



1 / 1 {i + 3wy 

-- < w < -, , , 

3- 3 8(3u;2 + 2ti; - r 



< /3i < 



(37) 



which means that, if Point D exists, it is always stable. The left panel of Fig. E] shows 
the region of parameters corresponding to Point D. We summarize the existence 

and stability of points C and D in Table |TT1 




Point C 



^ \ 

0.20 0.25 0.30 0.35 0.40 

w 



FIG. 3: Regions of stability in the {w^fi\) parameter space for Point D (left panel), and the 
evolution of Points C and D with the decreasing of w (right panel). In the right panel, A = 0.6 
and = —\ 
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FIG. 4: The phase transition from a stable state to an inflation by assuming a slowly deceasing 
equation of state(tt;(t) = 0.280 — O.OOlt) for the HL cosmology with the initial conditions a(0) = 
0.5 and d(0) = 0. The parameters are set as A = 0.6 and ^1 = — 1 



Thus, if the cosmic scale factor satisfies Eq. fl30|) initially, and w and /3i lie in the 
region given in ( l37l) . the Universe can stay at a stable state past-eternally and undergo 
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an infinite oscillation. If w evolves in such a way that w and Pi satisfy Eq. ( 13T|1 . then the 
stable critical point D coincides with the unstable one (Point C) and becomes unstable. 
As a result, the universe goes out of the stable state and enters an inflationary phase 
naturally. A particular case which realizes a phase transition from a stable state to an 
inflation era is shown in Fig. HI So the big rip singularity may be avoided successfully in 
this case. 

TABLE II: Summary of the critical points and their stabihty in the SVW HL cosmology with 





Existence 


Stabihty 


Point C 


w < < /3i < 1 

1 / ... / 1 {l+3wf / / 3 
3-^-3' 8{3w'^+2w-l) - Pi - 8 


unstable 


Point D 


1 ^ ,,, / 1 (l+3u>)2 ^ 5, / n 
3 < W < 3, 8[3w'-i+2w-l) < Pi < U 


stable 



2. /32^0 

Now, the Einstein static points satisfy Eq. (^7^, which is a cubic equation of a'^^. The 
solution of Eq. ( HT} is determined by the following expression 

A = - AAC , (38) 

where A = b'^-Sac, B = bc~9ad, and C = c^-3bd with a = 12(w-l)/32, b = 2{3w-l)Pi, 
c = —(1 + 3w) and d = {1 + w). 

• A > 0: there is only one real solution, which corresponds to only one critical point: 

Point E : -V = -^[b + n'^' + ^"2'^"] , (39) 

where Yi^2 = Ab + f{-B ± ^/A) 

• A < 0: there are three different real solutions. Thus, in this case, there are three 
critical points: 

Point F : = - — [b + 2VAcos{9/3)] , (40) 

Aa^g 3a 
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Point G : — ^ = -^[b- v^[cos(^/3) + ^3 sin(e/3)]) , (41) 



''Es 

Point H : = -^{h - VI[cos(0/3) - y3sin(0/3)]) , (42) 

where 6 = arccos (T) and T = ^^{2Ab - 3aB). 

• A = 0: Points G and H coincide since ^ = and thus there are two critical points 
(Points F and G with 6 = 0). 

At these critical points, the corresponding energy density p has the form 

with -4- given in Eq. f l39p -f l4ip . or fH21) . The conditions for these points to be physically 
meaningful are that p{aEs) > and a|;^ > 0. Since it is not an easy task to obtain analytic 
solutions to the existence conditions, we resort to numerical calculations and find that in 
order to satisfy the existence conditions, for Points F, G and H, /3i is required to be less 
than about |, while for Point E, there is no constraint on f3i. 

In order to show the regions of existence for Points E, F, G and H in the {w, (32) 
parameter space in detail, we chose —2, and 2 as three typical values for Pi. The results 
are shown graphically in Figs. OHl Figure shows the regions of existence for Point E 
with Pi = 2,0 and —2, Fig. [6l the region of existence for Point F with /3i = and —2, 
and Figs. [7] and [8] the regions of existence for Points G and H with /3i = and —2. 

In order to discuss the stability of these critical points, we need to calculate the eigen- 
values, which can be expressed as 

at each point. We find that the critical points F and H are always unstable, while critical 
points, E and G, are always stable as long as they exist. A summary of these critical 
points and their stability is shown in Table lllli Therefore, if the initial condition is such 
that the cosmic scale factor satisfies Eq. (!39|) or pil) . the big bang singularity can be 
avoided. To illustrate the stability of Points E and G visually, we plot the regions of 
stability in {w,/32) parameter space in Figs. [9] and [TOl respectively. 

14 



1 



-1 

A -2 
-3 
-4 
-5 



/J, =2 



2 4 6 
w 



^2 0.0 



8 10 





-1.5 



-1.0 -0.5 0.0 
w 



0.5 



FIG. 5: Regions of existence for Point E in the {w,l32) parameter space. The left, middle and 
right panels show the case with /3i = 2, and —2, respectively. 
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FIG. 6: Regions of existence for Point F in the (w,/32) parameter space. The left and right 
panels show the case with Pi = and —2, respectively. 
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FIG. 7: Regions of existence for Point G in the {w,/32) parameter space. The left and right 
panels show the case with /3i = and —2, respectively. 
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FIG. 8: Regions of existence for Point H in the {w,f32) parameter space. The left and right 
panels show the case with /3i = and —2, respectively. 



1 


-1 

^82 -2 
-3 
-4 
-5 



/J, =2 



1.0 
0.5 
P2 0.0 
-0.5 
-l.Ot 



P2 









\=-2 




r 







2 4 6 8 10 
w 



-1.5 -1.0 -0.5 0.0 0.5 
w 



-1.5 -1.0 -0.5 0.0 0.5 
w 



FIG. 9: Regions of stability for Point E in the (il',/92) parameter space. The left, middle and 
right panels show the case with = 2, and —2, respectively. 
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FIG. 10: Regions of stability for Point G in the {w,f32) parameter space. The left and right 
panels show the case with /3i = and —2, respectively. 
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Let us note that Points G and H move closer and closer as A changes from A < to 
0, and coincide once A = 0, which means that the stable Point G becomes an unstable 
one. Hence if the cosmic scale factor satisfies the condition given in Eq. (14T]) initially, 
the universe can evolve slowly from a stable region to an unstable one with the decrease 
of as shown in the Fig. [HI Because of the fact that this unstable state will lead the 
Universe to enter an inflationary phase, therefore, in this case, the big bang singularity 
can be avoided and a subsequent inflation can appear naturally. Note, however, that if 
the cosmic scale factor satisfies the condition given in Eq. ( 139|) initially, the universe can 
stay at the Einstein static state eternally and thus avoid the big bang singularity, but it 
cannot evolve to an infiationary phase with the evolution of w. 




FIG. 11: The evolutions of the stable Point G and unstable one H with the decreasing of w. 
TABLE III: Summary of the critical points and their stability in the SVW HL cosmology with 





Critical point 


Stability 


A > 


Point E 


stable if it exists 


A < 


Point F 


unstable 


Point G 


stable if it exists 


Point H 


unstable 


A = 


Point G or H with 6* = 


unstable 


Point F with 6^ = 


unstable 
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IV. CONCLUSION 



The Hofava-Lifshitz gravity is a power-counting renormalizable theory, which has an 
anisotropic scahng between space and time in the UV hmit, and thus breaks the Lorentz 
invariance. Applying this theory to cosmology, one finds that the Friedmann equation for 
a nonflat universe is modified by a ^ term. The SVW HL theory is a generalization of the 
original HL gravity by keeping the projectability condition but abandoning the detailed- 
balance one. This generalization introduces an extra ^ correction term to the Friedmann 
equation and modifies the coefficient for the term as compared with the HL theory. In 
the present paper, we study the infiuence of these correction terms on the Einstein static 
state. In the case of HL cosmology, if the cosmological constant A is negative, we find 
that there exists a stable Einstein static state. The Universe can stay at this stable state 
eternally and thus the big bang singularity can be avoided. However, in this case, the 
universe can not exit to an infiationary era. So the big bang singularity problem cannot 
be solved successfully. By making an analytical continuation of the model parameters, a 



positive A can be obtained 



18|. But, in this case, there is no stable Einstein static state. 



For the SVW HL cosmology, when P2 = 0, we find that there exists a stable critical 
point and an unstable one. If the cosmic scale satisfies the condition given in Eq.f l30l) 
initially, the universe can stay at the Einstein static state past eternally. With the decrease 
of w, the stable point and the saddle one move closer and closer. Once w reaches a critical 
value, this stable critical point coincides with the saddle one and there is only one critical 
point, which is unstable. Thus the Universe can go out of the stable state and then enter 
an infiationary era. Therefore, the big bang singularity can be avoided and a subsequent 
infiation can occur naturally. 

When P2 7^ 0, our results show that if the cosmic scale factor satisfies the condition 
given in Eq. pTl) and the equation of state w is larger than a critical value initially, 
the Universe can evolve from a stable region to an unstable one with the decrease of 
w. Therefore, in this case, the big bang singularity can also be avoided and an infiation 
can appear naturally. However, if the cosmic scale factor satisfies the condition given in 
Eq. (l39i) initially, although the big bang singularity can also be avoided, the Universe 



18 



cannot evolve to an inflationary phase with the evolution of w. 

Note added: While we are in the stage of revising the manuscript, the stability of 
Einstein static universe in a IR modifled HL gravity is discussed in 38 1. 
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